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1 dimensions as the candidates for the disoriented chiral condensate (DCC). 
The nonuniformly isospin-orientated solutions are shown to be related to the 
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1 Introduction 



The soft particle production in a very high energy hadron-hadron or nucleus-nucleus 
collision is an interesting phenomenon. Occasionally, the collision creates a large 
number of low energy (small p t ) particles, mainly pion quanta, initially populating 
in a small interaction volume and subsequently undergoing a rapid expansion. The 
perturbative QCD is not applicable in describing the dynamics since it involves a 
large number of quanta and the interactions are highly nonlinear. One may anticipate 
some novel dynamical feature of the nonperturbative QCD. Although there is some 
evidence that the p t distribution of these particles follows the scaling law with an 
effective temperature, it is not clear whether or not these low energy particles can 
actually thermalize so that their distribution can be described by thermodynamics. 
In fact, some deviation from the thermal distribution in the very small pt region, say, 
p t < 100 MeV, has been observed though data are poor in this region at present. 

On the other hand, it has been suggested, first by Horn and Silver that these 
low energy pions may be described by a classical theory. The number (N) of quanta 
involved is large, the quantum fluctuation is suppressed by 1/y/N. In addition, the 
low energy theorem on the pion-pion scatterings dictates that the quantum corrections 
to the scattering amplitudes are suppressed by a factor of p 2 /(47r/ 7r ) 2 . More recently, 
a scenario of disoriented chiral condensate 0, [§ suggests that these low energy pions 
may be out of equilibrium and undergo a quench following a chiral phase transition, 
and their interactions should be described by the classical chiral dynamics. 

In this paper, we shall determine the possible classical field evolutions that these 
low energy pions may follow based on the nonlinear cx-model. The advantage of 
the nonlinear cx-model over the linear a model is that the constraint of the vacuum 
expectation value on the fields is built in and that the pion fields always describe the 
massless modes irrespective of the vacuum orientation in the background. The a mass 
is taken to be infinity so the low energy structure of the theory is evident. We have 
obtained in an analytic form a class of classical solutions to the nonlinear cx-model in 
3 + 1 space-time dimensions as the candidates of the disoriented chiral condensate in 
QCD. Our general solutions have a transverse momentum distribution and need not 
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be subject to a boost-invariance constraint. The solution with a nonuniform isospin 
orientation is constructed by the chiral SU(2) x 577(2) rotation from a uniformly 
oriented solution. In the limit of a boost invariance and no transverse momentum, 
our solutions reduce to those of Blaizot and Krzywicki || . We study the distribution of 
the neutral pion fraction / for the pions that disintegrate from the disoriented vacua. 
We find that the distribution dP/df = l/(2y/f) holds for the uniformly oriented vacua 
and also for the boost-invariant vacua with an infinitely large uniform spread in the 
transverse direction, but it does not hold for the vacuum whose isospin orientation is 
nonuniform in space-time. However this distribution should be correct if one selects 
pions from within a small region in the 7/-kj_ plot event by event. This conclusion 
is reached through the analysis using the classical field theory method and also by 
studying the quantum pion states. 

We organize the paper as follows. In Sec. 2, we start with the analysis for the 
boost-invariant solutions with no transverse momentum. In Sec. 3, we make the ob- 
servation that all solutions with a nonuniform isospin orientation are obtained by the 
chiral rotations from a uniformly oriented solution whose energy is degenerate with 
the nonuniform ones. In Sec. 4, we give a general solution with a uniform isospin ori- 
entation, from which we can obtain the nonuniform solutions by the chiral rotations 
according to the prescription given in Sec. 3. The general solution has a nontrivial 
transverse momentum distribution and is not subject to the boost-invariance con- 
straint. In Sec. 5, we examine the charge distribution of the pions disintegrating from 
these disoriented vacua. The picture of classical field theory leading to the distri- 
bution dP/df = 1/(2a/J) does not apply to the nonuniformly oriented vacua except 
in the boost invariant limit with zero transverse momentum. For a general solution 
which has rapidity and transverse momentum dependence, the distribution holds only 
within each small segment in the y-k± plot. The charge distribution is also studied 
from the viewpoint of quantum multipion states, following Horn and Silver [[!]]. The 
modification of the distribution is attributed to the fact that for the general nonuni- 
form solution, more than one orbital state is available for pions to occupy so that 
there are many different ways to construct multipion states. 
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2 Boost-Invariant Solution in 1 + 1 Dimensions 



In high energy hadron or nucleus collisions, the configurations approximately invariant 
along the collision axis are of particular interest. We first focus on this class of 
solutions ignoring the transverse spatial dependence. We choose a nonlinear a-model 
as the dynamical model for QCD at low energy. The solutions that we obtain in this 
Section are equivalent to those of Blaizot and Krzywicki [|| though they are dervied 
in a slightly different way in order to clarify a relation between the uniformly oriented 
solutions and the nonuniformly oriented ones which plays an important role when we 
extend our argument to the more general case later. 

The phase and radial representation of the nonlinear cr-model is, 

E(x) = e iT - n ^ x \ (1) 

No matter what values the classical phase fields take, the state remains at the bottom 
of the potential valley because |S| = 1. This facilitates greatly the search for the 
DCC-type solutions which are realized at the bottom of potential well. Define the 
pion field 

tt(x) = f„n(x)0(x). (2) 

where n(x) is an unit isovector field obeying n(x) ■ n(x) = 1. Alternatively one defines 
7r(x) by a + it • tt = f n T, with the constraint o = J~~f% — ir 2 . In this case the pion 
fields are given by ir(x) = / T n(x) sin9(x). In either case, n determines the isospin 
orientation of the pion field. 
The lagrangian is given by 

£=S-tr(d^(x)d^(x)), (3) 

where £ transforms like E — > Ul^Ur under SU(2)l x SU{2)r rotations. In terms of 
6(x) and n(x), the lagrangian is 

£ = ^-(0^9 d»0 + sin 2 ed^n ■ d"n) + ^(n 2 - 1), (4) 

where A is a Lagrange multiplier. We will not include an explicit chiral symmetry 
breaking throughout this paper. The Euler-Lagrange equations are 

U0 = sin 9 cosflfyn- <9 M n, (5) 
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d^w^Od^n) = An. (6) 

The chiral SU(2)l x SU(2)r symmetry assures the conservation of the vector and 
axial-vector currents. In terms of 9 and n, the current conservation is written as 

<9 M (sin 2 #n x <9 M n) = 0, (7) 

<9 M (n<9^ + sin# cos#<9 M n) = 0. (8) 

The isospin current conservation (7) follows also from (6), while the axial- vector cur- 
rent conservation (8) can be derived from (5) by repeated use of (n • dn/dr) = 0. (4), 
(5), (6) and (7) are most general with no assumptions or approximations made. 

We consider a boost-invariant case in 1+1 dimensions where the fields 9(x) and 
n(x) are only functions of the variable r: 

r = Vt 2 - x 2 . (9) 

For a function only of r, a partial derivative <9 M /(r) is equal to (x^/r)df /dr. Fur- 
thermore, d^(f[r)d^g[r)) = (l/r 2 )(d(r 2 fg')/dr) where g' = dg/dr. The current 
conservation relations can be integrated into 

r sm z 9 n x — = a, (10) 

GST 

d0 . n n dn 
rn- — h r sm cos — = b, (11) 
dr dr 

where a and b are constant vectors in the isospin space whose magnitudes are denoted 
as a and b respectively. It is immediately obvious from (10) and (11) that a and b 
are orthogonal to each other: 

a±b. (12) 

The isovector field n(r) stays perpendicular to a as r varies. By multiplying (10) 
with n vectorially and using n • (dn/dr) = 0, one obtains 

dn _ a x n 

dr rsin 2 0' 1 j 
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a standard equation for a vector n to precess around a constant vector a. The pre- 
cession frequency |a|/rsin 2 # varies with the proper time r. The relations among a, 
b and n are illustrated in Figure 1. Squaring (11) gives 

('£) i+ - ws ( T $ i= *"- (14) 

Eliminating dn/dr from the these equations, one obtains the differential equation for 
»(t): 

where a = |a| and b = |b|. (13) and (15) combined contain the same information as 
the first integrals of the Euler-Lagrange equations for 6 and n so that we may proceed 
with the current conservation laws. 

(15) is analytically integrable into the most general boost-invariant solution for 
Q(t) in 1+1 dimension: 

cos#(r) = (6/k) cos (re In (t/t ) +0 o ), (16) 

where k 2 = a 2 + b 2 and cos^ — cos9(t ). Substituting 6{r) in the integral from 

the isovector current conservation, one obtains the general solution for n(r). Since 
n(r) and b both lie in the plane perpendicular to a, it is convenient to express the 
unit isovector field n(r) in terms of a single angle (5{t) measured from the direction 
of b: 

n(x)-b = 6cos/3(x). (17) 
The solution for cos f3(r) can be expressed in terms of 9(t) 

(18) 



cos (3{t) = 



K 2 



b\\a 2 sm 2 e(r)' 

The DCC configuration may be obtained from the boost-invariant solution multiplied 
by a step function 6(r 2 ): 

n^)Dcc = e^ x)T - nW 0(r 2 ), (19) 
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such that the causality condition is satisfied. Our solutions should only apply to the 
inside of the light-cone. Once 0(r 2 ) is inserted, there appears a source on the light 
cone at r = which triggers the formation of a DCC. The energy density S{x) of our 
solution is singular as we approach the light cone: 

S(x) = S. (^-)(a 2 + 6 2 ), (20) 

for both the uniform and nonuniform solutions. The isospin vectors a and b enter the 
energy density in the combination of a 2 + b 2 as required by SU (2) x SU(2) invariance. 
The lowest energy solution of a 2 + b 2 = is a trivial solution obtained from (16), (17) 
and (18) by taking the limit of a, b — > 0: 

9{x) = constant, n(x) = constant vector. (21) 

Blaizot and Krzywicki ||] expressed the pion fields by n = f n ii sin 9. If we choose our 
initial condition such that sin i9 = 0, that is, cos 6{tq) = b/n, our general solution 
given by (16) and (18) coincides with theirs. 

3 Chiral Rotation and Nonuniform Solution 

In this Section we focus on the relation between the space-time dependence of the 
isovector field n(x) and SU(2) x SU{2) invariance of the lagrangian. It is easy to see 
in (18) that our solutions have a uniform isospin orientation when the isospin vector 
a vanishes. When a = 0, (3(r) = (mod 27r), that is, n points to the direction of b 
for all uniform DCC's. Because of the SU (2) x SU (2) invariance, it is always possible 
to rotate the vector a by an appropriate axial rotation to the direction of the vector 
b. After the rotation the solution has a uniform isospin orientation and is degenerate 
in energy with the nonuniform solution prior to the rotation. It is not unfamiliar that 
if a system possesses some symmetry, a set of infinitely many new solutions may be 
obtained by making the symmetry transformation on a single solution. 
To be explicit in the present case, let us rotate a uniform solution 

E (ac) = e m{x)T - no , (22) 
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where n is space-time independent. Upon a general chiral rotation parametrized 
by Ul = e^ T ' nL and Ur = e lvT ' nR , the uniform solution is rotated into = 
Ul^o^Ur 1 . The transformed 9 and n fields are given by 

c B = (c£C v + (n L ■ n R )s ( s r ^jc + 

((no • n fl )c € s„ - (n x n L )s^c v + (n x n L ) • s € s„) s , (23) 

and 

ns e = (n L St;C v - n R c ( s v + (n L x n R )s ( s r ^c 

+ (^n Ci:C v + (no x n x ,)s € c 7? + (n x n R )c^s v 

+ ((n ■ n L )n R + (n ■ n^)n L - (n L • n fl )no)s{S„^s , (24) 

where cq and sq stand for cos 9 and sin#, respectively, and so forth, while c = cos# 
and s = sin6> . For an isospin rotation, we choose £ = r] and n L = n R . Then the 
rotated fields are 

cos#(t) = cos6 l o(r), 

n = (no • n L )n L - (n - (n • n L )n L ) cos 2£ + (n x n L ) sin 2f . (25) 

Since this is a global isospin rotation, the resulting field is another uniformly oriented 
solution with the same 9{x). For an axial rotation, £ = rj and = — n^, in particular, 
if n^, n R _L n , one obtains 

cos6 l (r) = cos2£ cos#o(t), 
n(r) sin0(r) = n sin# (T) + n L sin2£ cos9 (r). (26) 

The axial rotations turn a uniform solution into nonuniform solutions. We can actually 
show that our general nonuniform solution given by (16) and (18) is reproduced with 
a suitable choice of the rotation angle: 

tan2£ = a/6. (27) 
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Then (26) gives 



cos 6(t) 
cos (3(t) 



(P/k) cos 9 q (t), 

n ■ n = sin6 l o(r)/ sin#(r). 



(28) 



The result in (18) for cos/9 is obtained by solving the above equations. In this way 
we are able to obtain the nonuniform solutions from uniform ones by the axial- vector 
rotations. In the boost-invariant 1+1 dimensional case, we have obtained in Sec. 2 
all nonuniform solutions by solving explicitly the nonlinear differential equations. We 
have shown that they are all related by the axial rotations to the nonuniform solutions 
with the same energy density S(x). All possible solutions are exhausted in this way. 

4 General Solution Without Boost Invariance 

Even without boost invariance, it is straightforward to solve for the uniformly oriented 
solutions. Once we obtain a uniform solution, we can transform it into the nonuniform 
solutions by SU(2) x SU(2) rotations. The equations for 9(x) and n(a;) from (5) to 
(8) in Sec. 2 are also valid in the boost-noninvariant case. We look for the uniform 
solution in which n = constant so that d^n = 0. In this case the equation of motion 
reduces simply to 



The origin of space-time coordinates is identified with the collision point of the hadron 
collisions. The z-axis is chosen along the collision axis of the initial hadrons. Note 
that the meaning of variable r is a little different from the 1+1 dimensional case. The 
surface of r = lies outside the light cone with respect to the space-time origin except 
for the exactly forward and backward directions. With these space-time variables, 
□6* = becomes 



no = 0. 



(29) 



It is convenient to use the space-time variables 




(30) 




(31) 
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Since the differential equation is homogeneous, it can be solved by the method of 
separation of variables in the form 

9{x) = T(r)H(r])X( X± ) 9(r 2 - xi). (32) 

We solve for 9(x) inside the light cone, r 2 — xi > 0. 

For the transverse direction xj_, the general solutions are the Bessel and the Neu- 
mann functions. If we require that the solution be regular on the collision axis 
P — l x _i_| — 0, the Neumann functions are excluded. Note however that a singu- 
lar behavior means an infinite oscillation toward p — 0, not an indefinite increase, in 
terms of the pion field 7r = f w ii sin 9. One may also require that X(xl±) should not 
increase indefinitely as p — > oo. With these requirements, the parameter p 2 defined 
by Aj_X = —p 2 X must be positive. We choose p > 0. X(x±) is expressed in terms 
of the Bessel functions of integer order: 

oo 

X(x ± ) = CqJ q (pp) + J m {pp){C m cos m0 + D m smm(p), (33) 

m=l 

where Co, C m , and D m are the numerical coefficients to be determined by the boundary 
conditions. The magnitude of p determines a transverse size of a DCC and therefore 
a spread of the p t distribution of the final pions. Since a DCC will have an extended 
size in the transverse direction, the value of p is likely to be a fraction of f n or less. 
The spatial rapidity dependence H(rj) is simply solved 

H(r)) = cosh \i] or sinh Xrj. (34) 

The parameter A can be any complex number in general (If it is complex, one should 
take the real part of 9(x) at the very end). The special case A = leads to the boost- 
invariant solutions. For the approximately boost-invariant DCC configurations, the 
magnitude of A is much smaller than unity. The region of large values of r) corresponds 
to the forward and backward edges of DCC where energetic leading hadrons are 
moving outward while the small values of rj describe the cool central region. In this 
picture, it appears appropriate to choose H(i]) such that the pion density is higher at 
a larger rj than at a smaller 77. We therefore choose A to be real (and positive) in the 
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following. It should be emphasized however that our choice for a real A over purely 
imaginary or complex A is more for the convenience of the presentation. 
Given X(x.±) and H(r)), T(r) obeys the differential equation 



The solution is given by J\(/j>t) and/or N\((j,t). The main difference between J\(fir) 
and N\(/it) is their different behavior as fir approaches 0. In the limit of a boost- 
invariance, A — > 0, J\(fir) approaches unity while iV A (/rr) — > m(//r). If there is no 
transverse momentum, the latter approaches the uniform solution (a — 0, k — b) 
described in Sec. 2, while the former coincides with the lowest energy solution given 
in (21). Putting all together, one obtains the uniform solution in a complete form: 



where one may superpose these solutions in A, fi, and m. 

It is much harder to solve directly for nonuniform solutions when there is no boost- 
invariance constraint. Though some simple special solutions can be obtained by luck, 
finding all nonuniform solutions is a formidable task. In contrast, it is straightforward 
to perform SU (2) x SU (2) rotations on the uniform solutions. The rotation formulas 
(23) and (24) are most general and applicable to the boost-noninvariant case as well. 
Therefore the nonuniform solutions can be obtained by the chiral rotation from the 
uniform one in (36). 

An important question is whether we exhaust all nonuniform solutions by the 
axial rotations from the uniform solutions. In other words, are there any nonuniform 
solutions that cannot be rotated into a uniform one? If such a class of solutions exists, 
it would have some topological quantum number like a soliton. Note however that 
the solutions of our interest are time-dependent, and that their energies and actions 
are not necessarily finite. Unless these topologically nontrivial solutions exist, the 
uniform solution (36) and the SU (2) x SU (2) rotations on it exhaust all solutions. 




(35) 



9(x)n 



(aJ X (l2T)+bN X (/2T)) 




(36) 
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5 Pion Charge Distribution 



It has been predicted that the pions decaying from a DCC will show a distinct charge 
distribution when the charge ratio is plotted event by event in / = N n o / (N T o + N n ±). 
The distribution 

— - — (37) 
df ~ 2^7 (37) 

has been derived in two very different ways. The first derivation assumes that an 
isosinglet multipion state is created by the decay of a DCC Jl|. All pions decaying 
from a given DCC are assumed to occupy an identical orbital state that is determined 
by the spatial configuration of DCC. The Bose statistics allows only one isosinglet 
2iV-pion state: 

|2iVvr) = (2a{aL - aJaJ^O), (38) 

where a}± are the creation operators of the pions in the same single orbital state. 
Making a binomial expansion of the right-hand side at large N, one obtains a simple 
rule dP/df = l/{2y/J). It is later pointed out that the relative phase between a\_al_ 
and QjqQjQ IS inessential to the final prediction of dPj df || . 

The second derivation is based on a more intuitive picture in classical field theory. 
Assuming that the isospin orientation is uniform in space-time and that all isospin 
directions are equally probable, one obtains dP/dVt = 1/47T, where Q is the solid angle 
for an isospin direction in isospin space. Since the ir° fraction / is proportional to 
the square of the third component of the pion field, (/ oc cos 2 /3), one obtains again 



distribution (37). In this derivation, the interference effects are completely ignored. 

Let us examine whether or not this prediction remains valid for the nonuniform 
DCC's. In the first derivation, it is crucial that only one orbital state is available 
for pions and therefore the isosinglet state is unique: For two pions, the isosinglet is 
nothing but (2a\_al_ — aJa^O) by the Clebsch-Gordan coefficients. For four pions, 
the group theory alone would allow two isosinglets. One is to combine the from 
l n ® In — 02tt + 12tt with the other from l n ® 1^ = + 1%*- The other is to 
contract the \ 2 -k from l n <g> l n = + 12tt with the other l 2n from 1^(81^ = + 12tt- 
The Bose statistics forbids 1-2^ for two identical pions in the same orbital state so that 
only C^jrCSC^TrlO) = (2a) + a)_ — ajag) 2 ^) is allowed. This argument goes through for any 
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2N, leading to the |2iV7r) in (38). If there are more than one orbital states available, 
the four-pion singlet state would generally take the form 

|4tt) = (A(0 2n ® 2 ^ + B(l 2n <g> l 27r )) |0), (39) 

where the coefficients A and B are dynamics-dependent. There are increasingly many 
more isosinglets for 67r's, 87r's etc, as N goes up. In the above example, the A-type 
term and the 5-type term give quite different pion compositions: there is 7r°7r 7r 7r° 
in the A-type term, but no 7r°7r 7r 7r° in the 5-type term. In order to obtain the 
distribution dP/df = 1/(2^), there must be only the A-type term and nothing else 
in the 2Ntt state (N -> oo). 

One can construct explicitly the in state when the isovector field n(x) is nonuni- 
form in space-time. Let us parametrize the direction of n(:r) by the azimuthal and 
polar angles a(x) and (5{x) with respect to the isospin z-axis. To simplify our com- 
putation a little, we consider as an example a DCC whose isospin is nonuniform only 
in the polar direction f3, but not in the azimuthal direction by choosing a — 0. We 
shall use the representation n(x) = f w ii(x)sm6(x) instead of n(x) = f 7T n(x)9(x) for 
the following discussion since the former automatically incorporates the periodicity 
of E(x) in 9(x) — > 9(x) ± 27r. The Cartesian isospin components of the pion field are 

7Ti = f n sm9(x) sm(3(x), 

7T 2 = 0, 

7r 3 = f w sm9(x) cosP(x). (40) 

The DCC state is described by the quantum coherent state, up to an overall normal- 
ization 

\DCC(6,/3)) = exp(4( Se ^) + 4(s eC/3 ))|0), (41) 

where 

aldose) = J y^Mk^k^k, 

a\{s e c p ) = J V / 2^0 sc (k)4(k)rf 3 k, (42) 
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with tfi ss and <p sc being the three-dimensional Fourier transforms of f n sin 9 sin /3 and 
f n sin 9 cos (3 respectively. Unlike aj(k), the operators a[{soSp) and a\{secp) are not 
canonically normalized, but the normalization is irrelevant to the isospin structure. 
The | A^7r) projection of the DCC state is 

\Nn(9(x) P(x))) = l^JaKses^ + aUsec^yiO). (43) 

Under the assumption that the DCC's appear in the intermediate state with 1 = 
and the production processes conserve isospin, if one DCC can be produced, all other 
DCC's that are related to it by the isospin rotations can be produced with an equal 
probability. The isosinglet DCC state can be constructed from the state in (41) by 
integrating out the Euler angles over the entire isospin space. 

The Air state of an isosinglet DCC is obtained by averaging |47r(^(a;) P(x))) over 
isospin space. The computation is straightforward though a little tedious. Up to an 
overall normalization, the result is 



|4tt(7 = 0)) = ^(|at( SeC/3 )| 2 + \ a \s eS/3 )\y - 4|at( SeC/3 ) x at( SeS/3 )| Z J |0>, (44) 

where |a^| 2 = 2a^aL — ajao. The first and second terms in the right-hand side are 
the A- type terms in (39), while the last term is the B- type term. As we anticipate, 
the isosinglet An state of the nonuniform DCC is no longer of the form postulated 
in (38). For a uniform DCC, that is, f3(x) — > constant, ^(sgCp) and ^(sosp) are 
identical up to a factor (in the 1+1 boost-invariant solution in Sec. 2, f3(x) is so 
defined that f3(x) — > 0, namely a 1 " — > 0, in the uniform limit). Therefore, the 5-type 
term cannot exist for the uniform DCC's. Our construction of the isosinglet 4ir state 
and the existence of the B-type terms cast a serious doubt on the distribution for the 
nonuniform DCC's. 

Alternatively, let us study the problem by assuming that |DCC(6 l (x), (3(x))) with 
different 9(x) and (3(x) do not have the quantum interference with each other. It is in 
accordance with the classical field picture. For a large number of pions, ignoring the 
interference may be justified. The momentum spectrum of pion quanta (i = 1, 2, 3) 



2\ 2 
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decaying from a classical field is given J7[ 

= mm, ,45) 

where p~i(k, |k|) is the four-dimensional on-mass-shell Fourier transform of the pion 
source function Pi(x, t) defined by □7r(x) = p(x): 



Pi 



(k, |k|) = / pi(x, t) e-' ik x+i|k| 'dxrft. (46) 



It is convenient to perform the space-time integral using variables r, 77 and for 
which dx.dt = rdrdrjdiL^, and 

e =^ »=H^h)' k± ' (47) 

for the momentum variables. ( 45[) becomes 



d IV- r 2 

(2tt) 3 -— ft (r, ??, Xl ) e^^^^^^^^rdr^^X! . (48) 
dyd z k± J 

If a DCC is boost-invariant along the collision axis, p(r, 77, xj_) does not depend on 77. 
In this case, 77 is integrated out and the energy spectrum dNi/dy c? 2 kj_ is independent 
of the rapidity variable y, as it is well known. 

Let us look into the isospin structure of the Fourier transform of the source that 
enters the right-hand side of (48). After the space-time integration is performed, the 
isovector p is generally of the form 

p(k,|k|) = F(y,kJc(k), (49) 

where F(y, kj_) is an isoscalar, Lorentz-scalar function of k and of whatever parame- 
ters that characterize a DCC; e(k) is a unit vector in isospin space. The pion spectrum 
is simply 

For each fixed k, one may repeat the classical field derivation and reproduce 

dP 1 
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(51) 



However, it is clear that the pion number Ni no longer obeys distribution ([IT]) when the 
momentum k is integrated over. To illustrate this point, consider a toy DCC for which 
e(k) points to one direction for a half of the range of rapidity y and to another direction 
perpendicular to it for the other half of y. Such a DCC is not one of the solutions that 
we have obtained, but it serves to make a point. Since there is no way to align the two 
e(k)'s to the same direction by isospin rotation, there are no DCC's in this isospin 
family that emit only 7r°, even though all directions are equally probable in isospin 
space. For a family of nonuniform DCC's, dP/df is zero at / = (Centauro) and at 
/ = l(anti- Centauro) , and tends to bulge in the central region of /, unlike that for a 
family of uniform DCC's. Only if the uniform DCC's dominate over the nonuniform 
ones, can distribution (|3Tf) hold approximately. The abundance of the uniform DCC's 
has a measure zero relative to that of the nonuniform DCC's in the phase space of 
the rotation angles. Unless the production of the nonuniform DCC's by the initial 
hadrons is strongly suppressed for some dynamical reason, dP/df = l/(2\/J) cannot 
hold even approximately. The spectacular Centauro and anti-Centauro events will be 
far rarer than our naive expectation based on the uniform DCC's. However, there 
may be a chance to observe distribution (|5T|) by selecting pions of the same y and k^ 
within small uncertainties. 

A special case is a boost-invariant DCC in 1 + 1 dimensions where p(x) is r\- 
independent so that p(k) is y-independent. In this case, e(k) becomes a constant 
vector independent of k and distribution ([37j) follows even for the nonuniform DCC's. 
One the other hand, for the 47r state that we studied in this Section, we see nothing 
special about the boost-invariant nonuniform case with k^ = from the general 
boost-noninvariant case. Do two derivations contradict with each other? It is difficult 
to make a connection between the two arguments. In analyzing the \Ntt(I = 0)) 
state, the interference between different DCC's is essential while in the classical field 
analysis, each DCC state is not an eigenstate of isospin and the interference from 
different DCC's is completely discarded. Though the both methods have led to the 
same dP/df distribution, it is not clear how much similar or mutually compatible 
their physical pictures are. With this unsolved uncertainty, we state our conclusion in 
a less assertive way: if we follow the isospin analysis of |2iW), we see no mechanism 
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that leads to distribution ([H]) for the nonuniform DCC's. If we argue instead in the 
classical field picture, the distribution does not hold except for the boost-invariant 
DCC with zero k^. However, distribution (|37| ) should hold for pions which are selected 
from a small segment of rapidity y and transverse momentum k^. 
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Figure caption 

Figure 1: The direction of n(r) relative to a and b. As r varies, n(r) precesses 
in the plane perpendicular to the vector a. c = b x a. For the uniform solutions, n 
stays in the direction of b for all r. 
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